Scale-free brain functional networks 
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Functional magnetic resonance imaging (fMRI) is used to extract functional networks connecting 
correlated human brain sites. Analysis of the resulting networks in different tasks shows that: (a) 
the distribution of functional connections, and the probability of finding a link vs. distance are 
both scale-free, (b) the characteristic path length is small and comparable with those of equivalent 
random networks, and (c) the clustering coefficient is orders of magnitude larger than those of 
equivalent random networks. All these properties, typical of scale- free small world networks, reflect 
important functional information about brain states. 

PACS numbers: 87.18.Sn 87.19.La 89.75.Da 89.75.Hc 



Recent work has shown that disparate systems can 
be described as complex networks, that is assemblies of 
nodes and links with nontrivial topological properties, 
examples of which include technological, biological and 
social systems lj. The brain is inherently a dynamic 
system, in which the traffic between regions, during be- 
havior or even at rest, creates and reshapes continuously 
complex functional networks of correlated dynamics. An 
important goal in neuroscience is to understand these 
spatio-temporal patterns of brain activity. This Letter 
proposes a method to extract functional networks, as re- 
vealed by fMRI in humans, and analyze them in the con- 
text of the current understanding of complex networks 
(for reviews see 0,0:0). 

Figure 1 shows how underlying functional networks 
are exposed during any given task. In these exper- 
iments, at each time step (typically 400 spaced 2.5 
sec), magnetic resonance brain activity is measured in 
36 x 64 x 64 brain sites (so-called "voxels" of dimension 
3 x 3.475 x 3.475 mm 3 ). The activity of voxel x at time 
t is denoted as V(x,t). We define that two voxels are 
functionally connected if their temporal correlation ex- 
ceeds a positive pre-determined value r c , regardless of 
their anatomical connectivity [4|, |5| . Specifically, we cal- 
culate the linear correlation coefficient between any pair 
of voxels, x\ and x 2 -, as: 
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FIG. 1: (Color online) Methodology used to extract functional 
networks from the signals. The correlation matrix is calculated 
and then used to define the network among the highest correlated 
nodes. Top four images represent snapshots of activity and the 
three traces correspond to selected voxels from visual (VI), motor 
(Ml) and posterio-parietal (PP) cortices. 



r(xi,x 2 ) 



(V{x 1 ,t)V(x 2 ,t)) - (V(x 1 ,t))(V(x 2 ,t)) 
a(V(x 1 ))a(V(x 2 )) 



(1) 

where cr 2 (V(x)) = (V(x,t) 2 } - (V(x,t)) 2 , and (•) repre- 
sents temporal averages. 

Figure 2 shows the degree distributions of networks ex- 
tracted using this method. The data was collected while 
the subject was opposing fingers 1 and 2 during 10 sec- 
onds, and then resting during 10 sec. We find a skewed 
distribution of links with a tail approaching a distribution 
p(k) ~ &~ 7 , with 7 around 2. This power law is more ev- 



ident for networks constructed with higher thresholds r c 
(more correlated conditions). For decreasing r c , a maxi- 
mum appears which shifts to the right. Despite changes 
in parameters, networks remain clearly defined indicat- 
ing that the main conclusions are robust with respect to 
the selection of parameters. The small inset in Fig. 2 
shows the distribution of links of a network constructed 
from the randomly shuffled (in time) voxels' signal. This 
network displays a Gaussian degree distribution in which 
the mean and width depend on r c . The largest values of 
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FIG. 2: (Color online) Degree distribution for three values of the 
correlation threshold. The inset depicts the degree distribution for 
an equivalent randomly connected network. 



the correlation thresholds used to construct the random 
networks are usually extremely low (r c ~ 0.1) compared 
to that used to define the functional networks (r c ~ 0.7). 
Our data was also compared with values from a randomly 
re-wired network, where nodes keep their degree by per- 
muting links (i.e., the link connecting nodes z, j is per- 
muted with that connecting nodes fc, I) 6] (see below). In 
this control the degree of each node is maintained but all 
other correlations (including clustering) are destroyed. 

To test the generality of these findings the same analy- 
sis was performed in 7 subjects across 3 task conditions. 
During data acquisition 7] subjects perform on-off finger 
tapping with three different protocols. In one case they 
are instructed verbally to start and stop tapping, in the 
other one the start /stop cue is a small green/red dot in a 
video screen, and in the last one the start/stop cue is the 
entire screen turning green or red. The results are very 
robust across subjects and task conditions. In particu- 
lar, the average of degree distribution (see Fig. 3) shows 
a clear power law scaling decaying as p(k) ~ &~ 7 , with 
an exponent close to 2. Although a precise fitting is ar- 
guably difficult, we find that for r c =0.6 7 = 2, for r c =0.7 
is 2.1 and for r c =0.8 is 2.2. This power law, indicating 
that the functional networks are scale- free, implies that 
there is always a small but finite number of brain sites 
having broad "access" to most other brain regions. Those 
well connected nodes are comparatively much more nu- 
merous in these networks than in a randomly connected 
network. 

As shown in the bottom panel of Fig. 3 the average 
probability of finding a link between two nodes, separated 
at least by a distance A, also decays as a power law. The 
significance of the scaling with distance is unclear because 
of the well known extensive cortex folding, which makes 
linear distance a dubious parameter. 

The scale-free character remains unaltered even for 
tasks engaging different brain regions. This is already im- 
plicit in the aggregated data of Fig. 3 (top panel), but we 
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FIG. 3: (Color online) Average scaling taken from 22 networks 
extracted from seven subjects. Top Panel: Average degree distri- 
bution. The straight line illustrates a decay of k~ 2 . Bottom panel: 
Average probability of finding a link between two nodes separated 
by a distance larger than A (using r c = 0.6). 



further corroborated this feature by analyzing two radi- 
cally different brain states: listening to music and finger 
tapping. As shown in Fig. 4, although the topographic 
distribution of the functional networks is very different 
for the two tasks, they have similar scaling behavior. For 
comparison, the standard activation map derived with 
the generalized linear model 8] is also shown. 

Now we turn to describe statistical properties of these 
networks: path length and clustering. The path length 
(L) between two voxels is the minimum number of links 
necessary to connect both voxels. Clustering (C) is the 
fraction of connections between the topological neighbors 
of a voxel with respect to the maximum possible. If voxel 
i has degree fc^, then the maximum number of links be- 
tween the ki neighbors is ki(k{ — l)/2. Thus, if Ei is 
the number of links connecting the neighbors then the 
clustering of voxel z, Ci = 2Ei/ki(ki — 1). The average 
clustering of a network is given by C = 1/N d, where 
N is the number of voxels. Clustering was analyzed also 
with respect to degree. The average clustering over vox- 
els with the same degree C(k) = l/N k J2j={i\ki=k} Q> 
where the sum runs over the N k voxels with degree k. 



3 




Degree k 



FIG. 4: (Color online) Comparison for two tasks: Panels A and 
B correspond to a finger tapping task while C and D to listening 
to music analyzed with our method or the standard fMRI linear 
model. Colors in pictures of panels A and C code the number of 
links detected with our method, and those in panels B and D the 
activation map build with standard model yj- The link distribu- 
tions (lower panel) show that the networks for both tasks are scale 
free. 



Table 1 summarizes the results for the networks ana- 
lyzed showing the average values (n = 22 datasets) for 
each threshold (r c , first column) used to construct the 
networks. Listed are A/", C, L, the average degree (fc), 
and 7. The clustering (C ra nd) and path length (L ran d) 
values of an equivalent random network are also included 
for comparison. Note that as the threshold r c increases 
the total number of nodes N decreases substantially, re- 
sulting by definition in more correlated networks. As a 
result, the number of nodes with at least one link de- 
creases, and consequently the (k) value decreases as well. 
In all cases, the coefficient C remains four orders of mag- 
nitude larger than C ran d- Networks randomized using the 
rewiring described by Maslov et al. p also have cluster- 
ing significantly smaller than the raw data (the order of 
10 -2 ). This feature, together with the similarity of path 
length of the original nets and their randomized controls 
(L and L ran ^), is indicative of a small- world structure 
0, • This property is robust as it does not depend on 
parameter r c . 



r c N C L (k) 7 C r and L ran d 

0.6 31503 0.14 11.4 13.41 2.0 4.3xl0" 4 3.9 

0.7 17174 0.13 12.9 6.29 2.1 3.7xl0" 4 5.3 

0.8 4891 0.15 6. 4.12 2.2 8.9xl0" 4 6.0 



TABLE I: Average statistical properties of the brain func- 
tional networks. 
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282 


0.28 


2.65 


7.68 


NA 


0.025 


2.1 


Macaque VC 


32 


0.55 


1.77 


9.85 


NA 


0.318 


1.5 


Cat Cortex 


65 


0.54 


1.87 


17.48 


NA 


0.273 


1.4 



TABLE II: Previously reported statistics of relatively smaller 
networks. None of these networks is scale- free. 

To our knowledge, this is the first report on the topo- 
logical structure of a large scale brain network. Previ- 
ous studies employing these statistical analyses have been 
limited to the small data sets of C. Elegans Q, and two 
neuro-anatomical databases 

una, 

the macaque visual 
cortex [ll[ and the cat cortex 12] (see Table 2). These 
studies did not demonstrate scale-free features. Compar- 
ison with the previous two reports indicate the following: 
although clustering in the present study is smaller in ab- 
solute value, it is still orders of magnitude larger than 
the random case (10 _1 vs. 10 -4 ), while in the previ- 
ous reports the clustering of the experimental data was 
just one order of magnitude larger than the randomized 
controls in the best case. Interestingly, the average con- 
nectivity (k) in all cases is of the same order, despite 
the huge differences in networks' origins and sizes. Ac- 
cordingly, this consistency may reflect some constraint (s) 
inherent to network construction. These quantitative fea- 
tures show that the human brain network examined here 
has small world properties, a finding that was previously 
postulated [3, L|. 

Figure 5 illustrates the dependence of two important 
features upon a voxel's degree. The first is clustering, 
found in many cases to scale as C(k) ~ k~ a , an indica- 
tion of hierarchical organization PjJll4(. We see, instead, 
a relative independence of clustering from degree. The 
second feature is that a highly connected node tends to 
connect with other well connected nodes. As shown in 
the bottom panel of Figure 5, there is a positive cor- 
relation between the degrees of adjacent vertices. This 
correlation, also called assortative mixing, is not typical 
of biological networks, but rather is distinctive of social 
networks Transitivity in correlations contributes to 
increase artifact ually the clustering coefficient, using par- 
tial directed coherence or Granger causality [l^ in the 
future should clarify this. 

In summary, we report statistical measures showing 
that the functional correlations of the human brain form 
a scale-free network with small world properties and as- 
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FIG. 5: (Color online) Top Panel: Plot of clustering vs. degree. 
Bottom panel: Plot of a neighboring node's degree vs. degree illus- 
trates the assortative feature. Symbols represents individual data 
and continuous lines the average values for nodes with the same 
degree. (Same subject shown in Fig. 2, with r c =0.6). 

sortative mixing. While some of these properties have 
been informally discussed, this work is the first quanti- 
tative description of these large-scale topological prop- 
erties, as well as the first report of an assortative bio- 
logical network. The scaling laws demonstrated here are 
robust across parameters (Fig. 1), subjects (Fig. 3), 
and task conditions (Fig. 4), suggesting they are invari- 
ant properties of an underlying dynamical network. The 
present results complement the extensive work done in 
the context of brain functional and effective connectivity 
• The present approach has additional important 
implications. Namely, these studies can be extended to 
cases in which standard fMRI techniques cannot be used 
for lack of subject cooperation, (e.g., Alzheimer's pa- 
tients). Because scale- free complex networks are known 
to show resistance to failure, facility of synchronization, 
and fast signal processing 18], it would be important to 
see whether brain networks scaling properties are altered 
under various pathologies. In that regard, techniques for 
investigation of communities' structures should be 
useful to analyze these aspects. Work on models |2(]j . 
is needed to further clarify specific origins of the scaling 



laws. Overall, the network properties uncovered here, of- 
fer a novel window to investigate the dynamics of brain 
states particularly in cases of dysfunction. 
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